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Abstract 

We study analytically the general features of electroweak symmetry breaking 
in the context of the Minimal Supersymmetric Standard Model extended by one 
Higgs singlet. The exact analytical forms of the renormalization group evolutions 
of the Yukawa couplings and of the soft supersymmetry breaking parameters are 
derived to one-loop order. They allow on one hand controllable approximations in 
closed analytical form, and on the other a precise study of the behaviour of infrared 
quasi fixed point regimes which we carry out. Some of these regimes are shown to 
be phenomenologically inconsistent, leading to too small an effective /i-parameter. 
The remaining ones serve as a suitable benchmark to understand analytically some 
salient aspects, often noticed numerically in the literature, in relation to the elec- 
troweak symmetry breaking in this model. The study does not need any specific 
assumption on tan (3 or on boundary conditions for the soft supersymmetry break- 
ing parameters, thus allowing a general insight into the sensitivity of the low energy 
physics to high energy assumptions. 
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1 Introduction. 



Extensions of the Higgs sector of the standard model or of the minimal supersymmetric 
standard model (MSSM) ]I| , is a suitable framework to assess the phenomenology of the 
search for Higgs-like particles, and in a wider context, that of the supersymmetric partners 
of these particles. The next to MSSM, dubbed hereafter (M+1)SSM || g, g |, g, where 
one SU(3) C x SU(2) L x U(1) Y gauge singlet supermultiplet is added [0 to the MSSM, 
has attracted interest, initially as a framework for a natural solution to the so-called \i- 
problem ||] and later on as a source for interesting phenomenology which can differ from 
that of the MSSM. Due to the modification of the Higgs sector, the phenomenology and 
the upper bound on the mass of the lightest (observable) Higgs is altered ||, while the 
modified neutralino sector can lead to unconventional signatures for the sparticle searches 

Our main concern in this paper is the pattern of electroweak symmetry breaking 



(EWSB) and the dynamical generation of the [i parameter. As was initially noted in (lC 
and studied intensively flO|, PI, the v.e.v. of the gauge singlet field tends to be much larger 
than those of the two Higgs doublets, implying the tendency for the singlet chiral superfield 
to decouple from the other superfields of the theory. This means that, apart from some 
ranges of the parameter space where the singlino is the lightest supersymmetric particle, 
the features of the MSSM are basically shared by its minimal extension. However, the 
analysis of the issues in (M+1)SSM was mostly done numerically, or, when analytically, 
only for small tan/3 |TI| (no Yukawa couplings for the tau and bottom were considered). 
Furthermore this is often done in the framework of universality of the soft supersymmetry 
(susy) breaking parameters, or in some scenarios with mild non-universality which gave 
results comparable to the ones with universality ]12| . 



In view of the importance of such features (and possibly their generalization to an ex- 
tended singlet sector), it is important to attempt an understanding of the generic pattern 
for the dynamical determination of the v.e.v. of the singlet, without any specific assump- 
tions about the GUT-scale boundary conditions or the magnitude of tan f3. In this paper 
we will address this issue fully analytically and from two complementary sides: firstly 
the analytical evolution of the various Yukawa (and gauge) couplings and the soft susy 
breaking parameters, as well as the determination of the regimes of least sensitivity to 
initial conditions, i.e. in the vicinity of the infrared quasi fixed points^ IRQFP); secondly 
the study of the EWSB equations in compelling regimes and the interplay between the 
magnitude of the dynamically determined Higgs doublet mixing parameter jl, the exper- 
imental lower bounds on chargino masses, and the relaxation of universality of the soft 
parameters. 

The paper is organized as follows: In section 2 we recall the basic ingredients of the 
(M+1)SSM and introduce our notations. In section 3 we deal with the renormalization 
group equations (RGEs), give the analytical integrated forms of their solutions along 
lines similar to [O, [15| and classify the ensuing IRQFPs regimes. Four regimes are found 



3 another possibility of least sensitivity, not considered here, could be the occurrence of focus points 
at phenomenologically acceptable energy scales much like in the MSSM|13(]. 
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generalizing the MSSM case |TjJ. Some numerical illustrations of the IRQFPs regimes 
are given in section 4. Section 5 is devoted to the study of EWSB constraints and we 
conclude in section 6. More technical derivations and results are given in Appendices. 

From a different standpoint, it is worth keeping in mind potential difficulties that 
can arise in the (M+1)SSM in relation to the appearance of cosmologically problematic 
domain-wall solutions. We, however, bypass in this paper such problems and possible 
solutions to them |i~7| . Another interesting feature of the (M+1)SSM on which we will 
not dwell is the possibility to break spontaneously CP symmetry. In this paper we will 
will assume, without further reference, to be in regions of parameter space where CP is 
broken neither explicitly nor spontaneously IT3 . 



2 The (M+1)SSM. 

In this model, the Higgs sector is constituted by two Higgs doublets Hi and H 2 , and one 
singlet S. The spectrum, compared to MSSM is richer (one more CP even, CP odd Higgs 
field, and one more neutralino). Introducing as in the MSSM the matter fields (Q, T, B, 
E, L) the superpotential reads 

W = XSH X .H 2 + ^S 3 + y t TQ.H 2 + y b BQ.H x + y T EL.H x (2.1) 

where the dot product represents the SU(2) scalar product, and the superfields T (B) 
respectively, the left handed antitop(antibottom), E, the left-handed antitau, and Q (L) 
the left handed doublets for quarks (leptons). All the parameters in the superpotential are 
dimensionless, and a mass term is forbidden by a discrete Z% symmetry. This symmetry 
also prevents S to take large v.e.v. (< S >) and an effective /i parameter is generated 
(/} = A < S >) of the order of magnitude of 100 GeV. But, when S develops a v.e.v. this 
discrete symmetry is broken, and domain wall solutions appear. It is well known that 
such topological defects are excluded by cosmology. Possible solutions to this are given 
in [|17j. 

Supersymmetry breaking is parameterized by the so-called soft SUSY terms involving 
trilinear couplings (A : s), scalar and gaugino masses (m's and M's). In (|2.1|) Hi,H 2 etc. 
represent the superfields and now, although Hi, Hi etc represent its scalar component, 
Ai, A2 and A3 are the U(l)y, SU{2)l and SU(3) C gauginos respectively. 

£ soft = MiXiXi + M 2 X 2 X 2 + M 3 X 3 X 3 + (XA x SH 1 .H 2 + ^A K S 3 

+y t A t TQ.H 2 + y b A b BQ.Hi + y T A T EL.H x + h.c.) 
+m\\Hi\ 2 + m\\H 2 \ 2 + m 2 s \S\ 2 + m 2 Q \Q\ 2 

+m 2 T \T\ 2 + m 2 B \B\ 2 + m 2 E \E\ 2 + m 2 L \L\ 2 (2.2) 
Finally, let us write for later reference, the tree-level scalar potential in the neutral 
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Higgs sector including the F-terms from Eq. ( |2.1| ) , the D-terms and the contributions from 
the soft terms Eq.( [2.2|) , 

V = m\\Hl\ 2 + ml\Hl\ 2 + ml\S\ 2 + ^{\Hl\ 2 -\Hl\ 2 f 
+ \kS 2 + \H°H°\ 2 + \ 2 \S\ 2 (\H°\ 2 + \H° 2 \ 2 ) 

+(A X \SH°H° + A K ^S 3 + h.c.) (2.3) 

where g = \J (g 2 + ~g$)/2, gi,g2 being the gauge couplings associated respectively to 
U(1)y, SU{2) l and we consider only real valued couplings and mass parameters. 



3 Analytical solution for the renormalization group 
equation. 

3.1 Exact evolutions to one- loop order 

Using the notations a % = , i = 1, 2, 3; Y 5 = , j = t, b, r, and Y( X , K ) = where 
gi, g 2 and g 3 denote respectively the U(1)y, SU (2) l and SU (3) c gauge coupling constants, 
one can write down the one-loop RG equations as |3], |2Tf 



Mi = -biOiM, 

i I 

= - ( H c k iaiMi + a kiYiA^j 

i i 

t k = 2Y^c ki a i M i M i -J2akiY i (E i + A i A~) 

i i 

where k = t,b,r,X,K, . = d/dt, t=log Mq UT /Q 2 , the numerical coefficients a's, fe's and 
c's are given in Appendix A, and the E& are defined in Eq. ( |3.11| ). 

Though the RGEs for the Yukawa couplings Eqs.( |3.1| ) do not have explicit analytic 
solutions, they can be solved iteratively as it has been demonstrated in |[14]] through the 
use of some auxiliary functions. 

Together with the gauge couplings, the general solutions for the Yukawa couplings 



read 14 



a 



o 



Y, 



1 + hofit 
k 



1 + a kk Y k ° f* u k 
4 



(3.2) 



where the auxiliary functions u k are given by 



Uk{t) n fc (i + a^tiuj)^ (3 - 3) 

and 

E k (t) = Y[(l + b ia ° t t)- (3.4) 

8=1 

Specifying to the model under consideration one finds (see Appendix A for tabulation of 
the coefficients) 

_ /o r\ 

Ut ~ {l + QY°ju b )Hl + AY^u x )\ 

Ub = 1 ] r (3.6) 

(1 + 6Y t ° / u t ) e (1 + 4Y T ° / u T ) 4 (1 + 4Y A ° / u A ) * 

UT ~ (l + 6Y fe °/^(l+4Y A °/u A )i 

"a = 1 r— i 5" (3.8) 

(1 + 6^° / ut) > (1 + 6Y b ° J u b ) 3 (1 + 4F T / ur) * (1 + 6Y; / n K ) s 

m k = o- (3.9) 

(i+4r A n«A)t 

We also use the shorthand notation / to mean /q and drop out for simplicity any 
explicit reference to the scale t in all running quantities. Let us stress that ( |3.2|) give the 
exact solution to Y k and a, while the u^'s given in ( |3.5| )-( pl9|) , although solved formally in 
terms of the E^s and F fc °'s as continued integrated fractions, should, in practice, be solved 



iteratively. Yet, the important gain here is threefold: (i) as shown in |T3[], the convergence 
of the successive iterations to the exact solution can be fully controlled analytically in 
terms of the initial values of the Yukawa couplings. This means that one can in practice 
use truncated iterations of Eqs. (|3.5|) - (|3.9|) , say to first order, and thus obtain very good 
analytical approximations to the exact solutions (ii) the structure of the solutions is 
explicit enough to allow a thorough study of some limiting regimes as we will see in 
the next section (iii) furthermore, all these nice features will be naturally passed to the 
solutions for the soft SUSY breaking parameters since the latter will be obtained from 



through the method of JT5 



To obtain the solutions for the soft parameters one starts from those of the couplings of 
the super symmetric rigid theory, and makes the substitutions 



a -> a i (l+M i r) + M i fj + 2M i M i r ] fi), 

Y k -> Y k {l-A k r]-A k fi+(Z k + A k A k )rifi), (3.10) 
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where the Mi's are the gaugino masses, the A k s the scalar soft trilinear coupling constants 
and Sfc are the following combinations of the soft masses 



£t = rriQ + m\ + m|, E b = m\ + m% + ml, S T = m| + m| + m^, 

T,\ = m\ + wij + m\, S K = 3m| (3-H) 

and 77 = 8 2 and 77 = 9 2 are the spurion fields depending on the Grassmannian parameters 

e a ,h (« = 1,2). 

Performing the substitution (|3.10| ) in (|3.2| ) and identifying the coefficients of the re- 
sulting polynomial in 77 and fj, the linear term in 77 gives us the solution for Mj and A k 
and the 7]fj terms the ones for (For simplicity, we do not consider here CP-violating 
effects and take all the soft parameters to be real valued.) The resulting exact solutions 
look similar to those for the rigid couplings 



Mi 



1 + 



a „ 1 A k/ Y k ■+akkju k e k 

l/Y k u + a kk Ju k 
v _ , A 2 M9 A (Al) 2 /Y k °-i:l/Y k + a kk fu k Z k 

where the new auxiliary functions e k and have been introduced and are given in Ap- 
pendix B, Eqs.flETl]) - (p31) . 

In the particular case where Y\, = Y T = Yx = Y R = Eqs. (|3.2|) - (|3.12|) reduce to the 
exact well known solutions in the "small tan /3" regime. 

3.2 Large Yukawa regimes 

In this section, we study various regimes of large Yukawa couplings, and their incidence 
on the solutions written above. In ||16|| , iterated and analytical solutions to the RGE of 
the MSSM allowed an extensive study in the Infra-Red-Quasi-Fixed Point (IRQFP) fl9H 
regime, i.e., a regime where Yf =thr go to infinity. In the case of the (M+1)SSM, we want 
to see where and how strong could be the influence of the singlet on the MSSM solutions. 
For that, we have studied 4 different regimes including the two new Yukawa couplings 
Y® and Y K °, always considering Y i ( L tbT — > 00 in order to compare directly with the MSSM 



case If: : 



regime 1 : Y A ° and Y® are finite 

regime 2 : Y A ° is finite and Y K ° goes to infinity 

regime 3 : Y A ° goes to infinity and Y K ° is finite 
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• regime 4 : Y A ° and Y® go to infinity 

Let us make here some comments about the meaning of "infinite" initial conditions. 
Formally this means that we reach a Landau pole, and that the corresponding low energy 
values are at the edge of the triviality bounds. This clearly implies that perturbativity 
breaks down somewhere between the low and high (presumably GUT) scales. In practice 
though, we have checked numerically (see section 4) that we reach the effective fixed point 
(EFP) behaviour very quickly. A value of Y® of 0.1 is already in the EFP regime. So, 
there is no problem concerning the perturbativity at high scale. (In terms of the Yukawa 
couplings of the Lagrangian Eqs. (|2.1| ) - fl2.3|) we never take initial values larger than 5 
which corresponds to Y's < 0.16). 

When some of the Yukawa couplings become large, one is tempted to drop altogether 
the l's in the corresponding (1+Y® / itj)'s appearing in ( |3.2|) and (|3.5|)-(pT9|) and to expect 
a typical limiting behaviour for the w's of the formf] 

u°° ~ j—rr- (3.13) 

where Y° is a large Yukawa coupling. 

However, the situation is not so simple because of the implicit dependence of each 
on the full set of w's in a continued fraction like way, especially when not all the Yukawa 
couplings become simultaneously large, as is the case in some of the regimes we consider 
here. Although it is easy to understand intuitively the validity of (|3.13| ) if a similar form 



is obtained at some n th order iteration of the truncated approximation to Eqs.( |3.5| )-(|3.9|). 
and provided that the pi's verify < pi < 1 order by order, a more careful study is 
required to control the magnitudes of these powers. 

Indeed, in contrast to the MSSM case jnj, some of the p^s can be larger than one due 



to the singlet sector Eqs.(|3.8|), (|3.9|) . Technically, one has to solve ( |U.4[ ) of Appendix C to 
which the reader is referred for more technical discussions and details of the derivation. 

Here we give directly the final results for the running Yukawa couplings in the various 
IRQFP regimes: 



FP 



*Sb,r = 7 fp re S imes !> 2 > 3 and 4 ( 3 - 14 ) 

a H J U t,b,r 
yO M FP 

(yU)p A ~ regimes 1, 2 and 4 
Y™={ (3.15) 

regime 3 



4 /'< 



FP 



,FP 



4 note also that this approximation does not work at scales too close to the initial (GUT) scale because 
of the term L Ui, but it does hold for a typical electroweak scale t ~ 66. 
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FP 



Y° 



1_ 

61 



regime 1 



regimes 2 and 4 



regime 3 



(3.16) 



where the uf F are related to w°° through 



uT 



u 



FP 



(3.17) 



and depend, as well as the p^s, on the regime under consideration (see Appendix C). 
Let us stress several points here : 

(i) the solutions for the Y^F tbT have the same form as in the MSSM fL6f . Nevertheless, 
the effect of the singlet is implicit in the recursive solutions for the uf F Eqs. (|C.6|) - ( |C.10| ). 

(ii) Y A FP ~ in the regimes 1, 2 and 4. This will be important for the electroweak 
symmetry breaking discussed later on. 

(iii) We have an exact analytical solution for Y K in the regimes 1, 2 and 4, as a function 
of its initial value at the GUT scale, Y" K ° (see the numerical analysis for a more detailed 
discussion). 

(iv) It is also important to emphasize that the w FP Eqs. (|C.6|) - ( |C.10|) depend only 
on the ratios of the large initial values of the Yukawa couplings. However, even this 
dependence drops out completely in Eqs. (|3.14j) - ( |3.16| ), so that the initial conditions are 
completely screened in the IRQFP regimes as expected. Only a sensitivity to the initial 
Yukawa couplings that are not large may still occur, like in regime 1 for Y K FP , Eq.( |3.16| ). 

(v) A comment is in order here about the difference between the IR quasi fixed points 
we discuss and the exact fixed points studied in ||. For one thing, the latter exact fixed 
points, actually exact fixed ratios, exist only in reduced couplings configurations where 
all gauge couplings and all Yukawa couplings but the top are neglected, while the quasi- 
fixed points we study are valid without this approximation. Furthermore, the IRQFP's 
are more likely to influence the evolution from the GUT scale to the electroweak scale 
than are the exact fixed points Let us note, however, that in one or the other of 
our four IRQFP regimes we find either k or A to be vanishingly small, Eqs.( ^.l5| ), ( |3.16| ), 
similarly to the case of two among the three exact fixed point regimes determined in ||. 
Nonetheless, the latter two regimes were found to be infrared repulsive ||, while as one 
can infer from the structure of the denominator in Eq. (|3.2j ), the top down evolution of 
the Yukawa couplings tends generically always towards the IRQFP's behaviour. 



To find the IRQFP behaviour of the soft parameters Ai, Ej one can either perform the 
substitutions Eq.(|3TTq) in Eqs.([TI§ - ( j37g ), ( CT - QCTTCD , Q, or study the large 
initial Yukawa limit directly from the general solutions Eq. (|3.12|) . Denoting the auxiliary 
functions of Eqs. flB.ll) - ( |B.5| ) by e°° in this limit, we obtain for the A's 
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A™ 



00 + 



oo 



,FP 



A: 



= p in all regimes 
regimes 1, 2 and 4 

regime 3 



FP.FP 
A e A 



(3.18) 



(3.19) 



-e FP + k 



-e2° + 4° 



where 



regime 1 

regimes 2 and 4 
regime 3 

FP„FP 



(3.20) 



FP 



FP 



(3.21) 



and the e°° are defined in Appendix D Eqs.( p.2|) - ( |D.6|) . These auxiliary functions de- 
pend explicitly on the initial values A Q k . In some IRQFP regimes, this dependence cancels 
out in the running parameters A°°. Such an independence occurs in all the regimes for 
A tt b, T Eq. (|3.18| ) as was the case in the MSSM |U|], in the regime 3 for A\ Eq. Q3.19Q and in 
regimes 2, 3 for A K Eq.( |3.20| ). In such cases we have re-expressed the results in terms of 
new quantities ef p which are independent of the initial conditions. The dependence on 
initial conditions is explicited further in Eqs.( p.8|) , ( p.9| ). 

Similarly, denoting by £f° the auxiliary functions £j Eqs. ( |B.7| ) - ( |B.11| ) in the large Yukawa 
limits we obtain for the S's 



yoo 



^FP 



sf P + (AD 2 + 2er A 



FPtFP 



r bfcl 

FP\2 , r. FP A FP J U i Si 



FP 



yioo 
A 



^oo + ( A ooy + 2A^e^ + S° - (A° x ) 2 regime 1, 2 and 4 

P FP cFP 

e A FP + (A?)* + 2e£MP - regime 3 



(3.22) 



(3.23) 



FP 
A 



yO„/4()\2_ fiV f too 

C + {A™f + 2A-e- + C + ~ ( V J regime 1 



go + {A ooy + 2A o* e oc + E o _ {A oy 



regimes 2 and 4 
regime 3 



(3.24) 
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where again £f p have been used instead of £f° for those cases where the dependence on 
initial conditions is absent in the running E's. The reader is referred to Appendix D for 
an extended discussion of the relation between ef p ,£f p and e°°,£°° and for the explicit 
dependence on initial conditions Eqs. flDJj) , QD.19Q Some remarks are in order: 



(i) A\ depends on the initial condition A® in the regimes 1, 2 and 4, but in these 
regimes Y\ ~ Eq. (ft.lS ). Similarly, A K depend on the initial values of A® in the regime 



3, but in this regime Y K ~ (|3.16|) . So, we can conclude that, in every regime, the 
running of the combination Air^jYAi, i — X,k, the one present in the Lagrangian, is 
indeed screened from the GUT-scale initial conditions. 

(ii) We have the exact analytical dependence of the soft terms on the initial conditions 
A°,E? (c.f. Appendix D). 

(iii) Whereas the dependence of the £°° on the products A® A® and E° is rather com- 
plicated, the E°°'s depend only on the E°, except for regime 1 where A° K is present, see 
Eqs.( p.l8| ), ( p,19[ ). Moreover, this dependence is exactly the same as the dependence of 
the A°° on the A -. 

L J 

(iv) It is worth noting that the abovementioned sensitivity of to Y" K °, A^ in the 
regime 1 disappears in all the soft masses, except for the singlet soft mass (ms), see 
Eqs. flB.13l ) -( |B.20 ), ( D.36| ). Moreover, in regimes 1, 2 and 4, where the soft mass of the 



singlet separates from all the others, Eq. ( |D.35| ), the dependence on initial conditions for 



the soft masses is exactly the same as the one found in Jl6| in the case of the MSSM. 
Finally, let us note that it does not seem possible to give the explicit dependence of the 
A and the E's on gauginos mass initial conditions because the latter come always in scale 
dependent contributions, Eqs.( |B.6p , ( |B.12| ). 

A further point should be made here about the evolution of the soft scalar masses since 
we address the most general situation beyond universality. In fact, to solve exactly the 
RGE for the soft masses, we have to consider the complete equations, including a U(l) 
induced "trace term" 5*, 

K) = /.(Ml M 2 , M 3 , Ei, AA + T iai S (3.25) 
where /i(M li2j3 , Ej, Ai) are defined as 



ft = -a 1 M 2 1 +-a 3 M%-2Y t (Z t + A 2 t ) 

4 16 
f b = -a 1 M 2 + —a 3 Mi-2Y b (Z b + A 2 b ) 

1 16 
f Q = -a l Ml + Za 2 Ml + —azMl-Y t (Y> t + A 2 )-Y b (£ h + A 2 b ) 

f E = 4a 1 M 2 -2Y T (Z T + A 2 T ) 

f l = a 1 M 2 + 3a 2 M 2 -Y T (J: T + A 2 T ) 

f Hl = ai M 2 + 3a 2 M 2 -Y x (J: x + A 2 x )-Y T ^ T + A 2 T )-3Y b ^ b + A 2 b ) 

f Ha = aiM 1 2 + 3a 2 M 2 2 -r A (E A + ^)-3r t (Ei + ^) 
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f a = -2Y A (E A + ^)-2Y l6 (E K + ^) (3.26) 



S = {m 2 Q — 2m\j + m 2 D — m\ + m 2 E ) + m\ — m\ (3.27) 

generations 

and 

T% = {1/3, -4/3, 2/3, -1, 2, -1, 1, 0} (3.28) 
with i = {Q, U, D, L, E, H h H 2 , S} respectively. From Eqs. (^25|) , (|3~26|) , (|3~27| ) and ( ggp 



one sees that S 1 oc aiS*, so that if S 1 vanishes at some scale to it will vanish identically at 
any scale t. For instance S = at any scale in the case of universality of all soft scalar 
masses. More generally, S can still be ignored even when universality is relaxed provided 
that the initial conditions are such that S(t ) = 0. This simplifying configuration was 
taken up in \W\. Generically, however, one should solve Eq.( |j.25| ) keeping the trace term. 



This can be easily done by writing the solution for the soft masses as 

m] = {m 2 ) f + (m 2 ) Tr (3.29) 

where (m 2 )f is solution of the equation (m?) j = fi(M\, M 2 , M 3 ,£j, Ai), and ijn 2 ) Tr = 
TidiS. The (m 2 )f are then given by Eqs.( p.l3|) , ([B.20| ), and {m 2 )T T is just the solution 
of a linear differential equation system. We obtain 

(m 2 ) Tr = US {{1 + foaa?t)i - 1) (3.30) 

with 



^ = ^{1,-4,2,-3,6,-3,3,0} (3.31) 



where S = S(t = 0). 



4 Numerical analysis 

As we stressed before, one of the advantages of our approach is that we can extract the 
exact sensitivity to initial conditions in the quasi-fixed point regimes. In this section, we 
study numerically the solutions for the Yukawa couplings Y, the trilinear couplings Ai, 
and the soft parameters m,, in the four IRQFP regimes we have identified in section 3. We 
compare our analytical solutions to the results obtained by a purely numerical resolution 
of the RGE's, relying on a FORTRAN code which evolves the relevant parameters from 
the GUT scale to the electroweak scale through an algorithm similar to the one used in 

Ell- 
in the case of the Yukawa couplings we will restrain ourselves, for the sake of illustration, 

to the study of Y K . In Fig. la we plot Y K at the EW scale (t ~ 66) as a function of its 

value at the GUT scale, for different values of Y® and for Y t ° bT = 0.1. Since Y K ° is varied 

only in a range of small values, we reach in the same graph regime 1 (finite Y A ° and Y K °) as 
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well as regime 3 (infinite Y® and finite We have also plotted in the same figure the 

analytical solution (|3.16| ). We thus see, in the small Y® region, the numerical agreement 
with the the expected behaviour Y K (t = 66) ~ 1+6 y 0f of regime 1, and Y K (66) ~ for 
large values of Y" A ° in accord with regime 3. Fig. lb shows the variation of Y K (t = 66) 
for bigger values of Y® (regime 2 and 4). We clearly see that Y K tends to the analytical 
solution -^— t (~ 0.0025 for t = 66) for small values (regime 2) or high values (regime 4) 
of Y®. Evaluating the uf p, s up to the third order iteration with a procedure similar to 
that in fl6|], we also found very good numerical agreement between fl3.14|) , (|3.15| ) and the 
FORTRAN code output. 



We note finally that our results for Y K in regimes 2 and 4 are in perfect agreement with 
the numerical illustrations of || , while Y\ was found to differ slightly. The reason for this 
numerical difference can be traced back to the fact that Y K is completely decoupled from 
the t, b, t sector, while Y\ is closely tied to it. Whence an important influence of the Y t ^ T 
on the running of Y\, which modifies significantly its low scale value. 
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Fig. l.a Fig. l.b 

In figures 2a, 2b we plot the coefficients of the various initial conditions A®, (i = t, b, r, A, k) 
at the GUT scale which enter the running A\ at the electroweak scale (t = 66) as a function 
of Y x °. Since Y® is varied in a large range of values and the illustrations are given for one 
very small and one very large value of Y" K °, Fig. 2a, 2b cover all four regimes. Furthermore 

5 As far as the numerics are concerned, finite Y® means Y® ~ 10~ 2 and infinite Y® means Y® ~ 10 _1 . 
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the initial conditions for YfL t b T are fixed to the (large) common value 0.1. Thus when Y® 
runs from to 5 one observes the transition from regime 1 to regime 3 in Fig. 2a and from 
regime 2 to regime 4 in Fig. 2b. [We should keep in mind that the very large values taken 
by Y" A ° or Y® are only for the sake of numerical comparison in the deep IRQFP regimes at 
the electroweak scale. At much higher scales they eventually lead to perturbatively non 
reliable results.] The trend of the plots is in perfect agreement with what is anticipated 
from Eq.dETl). In Fig. 2a one retrieves the coefficients of regime 1 in the region of small 
Y" A °. For large Y" A ° all the coefficients vanish asymptotically as expected in regime 3. [The 
asterisks (*) at the right of the graph represent the values of the coefficients for y A ° = 10.] 
In Fig 2b where Y® is taken very large (= 20) we observe the evolution as a function of 
Y® from the regime 2, to the regime 4. The weaker sensitivity to Y® (as compared to Fig. 
2a) corroborates here the fact that the two regimes 2, 4 lead to identical coefficients in 
this case, see Eq.( 
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We have made the same analysis for the soft terms £, and m ; and obtained numerical 
results in total agreement with our analytical expression (Appendix D). We just present 
here in Fig. 3a, 3b, respectively the dependence of (?t?.q) 2 and (mi) 2 on the initial GUT 
scale values of all soft squared masses (not including the trace contribution of Eq. ( |3.30| ) , 
for a small value of Y" K ° (0.001), and for Y® running from to 5. One sees the transition 
from regime 1 to regime 3 with increasing Y®. The dependence on the soft mass (mg) 2 
vanishes towards regime 1. 
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Finally we illustrate in Fig.4 the behaviour of Y\, Y K in all four regimes, to stress the 
fact that these regimes do set in effectively, well before that the initial conditions Y®, Y® 
become infinitely large. 
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Fig. 4: Y x (full lines), Y K (dotted lines) 



5 Electroweak symmetry breaking 

In this section we consider the impact of the four IRQFP regimes on more phenomeno- 
logical issues. It turns out, as we will show hereafter, that regimes 2 and 4 are consistent 
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with the requirement of EWSB only for very small \fi\. These regimes are thus disfavored 
or excluded from present experimental exclusion lower bounds on the lightest chargino 



mass E3, 23]. Regimes 1 and 3 do not suffer from such features and correspond to viable 



configurations of least sensitivity to initial conditions. 

The study is carried out at the tree-level. [Some of the conclusions will remain true if 
one-loop corrections are included, but we will not discuss the issue further in the present 
paper.] The EWSB conditions involving the three Higgs v.e.v.'s are obtained from Eq. ( j2.3|) 
in the form: 

h x [ml + X 2 (h 2 2 + s 2 ) + \f(h\ - h 2 2 ) + (A\\s + nXs 2 )^-] = (5.1) 

2 il\ 

h 2 [ml + X 2 (hj + s 2 ) + \g\h\ - h\) + (A\\s + nXs 2 )^-] = (5.2) 

2 h 2 

s [m 2 s + X 2 (h\ + h\) + 2k 2 s 2 + A k ks + (— + 2«)A/i 1 /i 2 ] = (5.3) 

s 

where 





< H x >= "M < H 2 >= " < S s (5.4) 



Here hi, h 2 and s are chosen to be real- valued and positive, (see for instance 0] for a 
discussion of the freedom in the choice of the various parameters). 

In addition to Eq. ( |5.1| - 15.3| ) the requirement of correct electroweak scale reads 

,9 ,9 Ml / \ 

k l + k 2 = ~=f ( 5 - 5 ) 

From equations ( p.l| ), ( |5.2| ) and ( p.5| ), and provided that non of the v.e.v.'s is vanishing, 
one easily gets 



M 2 |M| + m|+AV 
1 g 2 M 2 +m 2 + m 2 + 2X 2 s 2 



M 2 \M 2 +m 2 + X 2 s 2 
2 g 2 M 2 +ml + m 2 2 + 2X 2 s 2 1 ' 

hh Ml (A\Xs + kXs 2 ) 

hih 2 = — —~ 772 (5.8) 

9 2 m f +m 2 + A 2^ +2s 2) 

Note that the above equations yield the familiar EWSB conditions of the MSSM 
when one goes to the limit A, k — ► 0, s — > oo with As, ks finite, with the identifications 
p, = As, B = A\ + ks, tan/5 = h 2 jh\, h\h 2 = M| sin2/5/^ 2 (of course with the proviso of 
Eq. Q5.3D which correlates dynamically s, thus jl, to the other parameters). 
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From Eq.( |5.8| ), one obtains^ 



with 



Y = a + b\ 2 (5.9) 
A 



A x 1 



1 M\ S 2(3 



2 P gi 



{c w M 2 z + 2m{) (5.11) 

(5.12) 



Equation ( |5.9| ) shows a linear correlation between k and A (in the regime of small A(<C g 2 )) 
even at low energy, the coefficient a depending on y t . This was noticed numerically in 
[p0|| for A and k° in constrained cases (universality, y t ^> y^...). 

On the other hand, the scalar potential at the electroweak symmetry breaking mini- 
mum is obtained from Eqs.( |2.3| ), ( (5.1H5.3Q and reads 



Vmin = -(KS + \h 1 h 2 ) 2 -\ 2 S 2 {h 2 l + h 2 2 ) 

-A x Xsh x h 2 - -A k ks 3 - - A M 2 z (h\ + h\)c w (5.13) 

Eqs.( |5.6| - 15.13| ) will be very useful when discussing the implication of IRQFP regimes. 

As previously stressed, although the IRQFP regimes correspond, strictly speaking, to 
some or all of the Yukawa couplings taking infinite values at some high energy scale, in 
practice we stay away from these unphysical configurations. Still, the typical IRQFP is 
essentially preserved as can be seen from Fig. 4, and gives a very good idea of the sensitivity 
to initial conditions. 

Bearing this in mind, we thus refer in the following to regime 1 as implying k,X g, 
to regimes 2,4 as implying A k, g and to regime 3 as k A, g at the electroweak scale, 
rather than the more strict behavior given in Eqs.( (|3.15|) , ( 3.16Q ). Since the discussion 



may get a little too involved, especially in regime 3, the reader interested only in the 
conclusions can go directly to the summaries in 5.2 and 5.4. 

5.1 Regimes 1, 2, 4: 



Note first that Eqs.(|5.6|, |5.7| ) do not have a solution if A = with s finite, unless m\ + m 2 . 



is fine-tuned to zero in which case there is an infinitely degenerate set (a valley) of solu- 
tions for hi, h 2 . This degeneracy is not lifted by one-loop corrections as can be seen for 
instance from the form of the top/stop contributions (jTOj and should thus be discarded as 



non physical. We now take A non vanishing but small enough to allow a reliable expansion 



we use from now on the shorthand notation, S2/3, C2/3 for tan (3, sin 2/3, cos 2/3. 
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to first order. Moreover, we consider separately the cases (a) s < m, A and (b) s ^> m, A 
where m, A denote generically the magnitudes of the soft masses and couplings. We will 
also denote by M so f t a generic value of the soft masses and couplings, or the smallest of 
these values. 



(a) s < m, A: 



In this case |A|s <C m, A and an expansion in small parameter is reliable. Performing 
this expansion in Eqs.( |5.6| - 15.8|) to second order in fi/M so f t (and first order in e) and 
eliminating hi, hi, the equation determining dynamically /t(= As) takes the simple 
form: 



^2, A X V(|Ml + m?)(|Ml + mI) 
^ + 2 , rn 2 ( 1 ~ e )^+ 1 F772 , —2 , =51 + °(( 



2 " Z "' Sn2 "' Z 1 - » ^ 3 ) = (5.14) 



M. 



so/t 



where 



A 2 M\ 
g 2 m\ + m\ 



+ 



(m 2 — m 2 ,) 2 



(M| + 2m 2 )(Mf + 2m|)|Mf + m 2 + m||(M| + m 2 + m|) 

1 K (i- £ ) 



^m 2 + m 2 ) A 



(5.15) 



(5.16) 



regimes 2, 4: 

Here |A| g. Solving Eq. (|5.14|) for fi leads in this case to 



fi = ±W--(mf + ml) 



(|M 2 +m 2 )(|M 2 +m 2 ; 
\| |M| + m 2 + m|| 



(5.17) 



which, apart from the consistency conditions (|M| + m 2 )(|M| + m 2 ,) > and 
A«(mf + t?t|) < 0, shows that tends to be generically very small being suppressed 
by the size of 



I A/ft I . Thus, even if one chooses the soft parameters sufficiently larger 
than the electroweak scale, so that the condition |A|s m, A dictated by regimes 
1,2,4 does not imply a priori small in comparison to the electroweak scale, then 
the dynamics of regimes 2, 4 will still lead to a vanishing fi. A vanishing fi implies 
a lightest chargino to be lighter than M\y and drops even much lower, for tan f3 > 
1, see Eq. (|5.38D. Such a configuration would be excluded, or at most marginally 



acceptable, by the LEPII lower bounds on the lightest chargino [22| |23], were it not 
for the fact that, since A is small in the regimes under consideration, unconventional 
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signatures due to displaced vertices can emerge in the case of the (M+1)SSM thus 
requiring a reanalysis of the experimental data ||. Even so, the theoretical upper 



bound ( 5.38 ) with vanishing /} will start conflicting with the conservative LEPI 
kinematical limit of M^/2 as soon as tan/3 > 2.27 or so. 

regime 1: 

This regime has a crucial difference with 2 and 4, namely that k and A can be 
generally of the same order. In this case the behavior given in Eq. (15.171 ) is no 
more valid. One can of course still solve readily Eq. (|5.14|) keeping in mind that 
<C m, A in the regime under consideration. Since the fx independent term in 
Eq.( [5.14D is generically of order 1/2 it is straightforward to see, taking into account 



Eqs. (|5.15| , |5.16|) that a consistent fx requires that A\ be very large compared to 



\Jrrii + m^, at the electroweak scale. This suggests that, generically (i.e discarding 
fine-tuned cancellations in m\ + m^), regime 1 disfavors solutions with s < m,A 
if the relevant soft parameters are of the same order of magnitude at some initial 
scale. 

(b) s > m, A: 

In this case | A|s > m, A. 
regimes 2,4: 



Eq.(|5.8|) takes the form 



k M 2 1 

^-A*¥ (^fTW + 0(AO) (5 ' 18) 

that is, h\h2 can be made arbitrarily large in the deep 2,4 regions ( ~ — > ). This 
is however in contradiction with the behavior dictated by Eqs. ( |5.6| , |5.7| ) in the same 
regions and leads to an inconsistency, even if M§ + m\ + m\ + 2fx 2 is artificially 
fine-tuned to 0(A). 

regime 1: 

Here no simple expressions can be derived and one would have to solve the full- 



fledged higher order polynomial equation for fx combined from Eqs.fl5T£- 



5.2 Summary for regimes 1, 2, 4 

We have shown that the IRQFP regimes 2,4 can be consistent with EWSB only when s 
is of the order of the soft masses, in which case \(x\ becomes too small to be consistent 
with present limits on chargino masses (or at best marginally consistent if tan/3 < 2.3 
when only conservative LEPI limits are considered). On the other hand we found that for 
dynamical reasons, and independently of any phenomenological considerations, regime 1 
can be consistent only for large s (~^> m, A) . This last point is one ingredient for the 



explanation of the numerically established large values of singlet Higgs v.e.v., |TOfl. We 
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will come back later to these features. Let us also note that the above conclusions were 
drawn without taking into account Eq. (|5.3| ). This equation can be viewed in this context 
as nearly correlating the two extra soft parameters m s , A K and enters the game as a further 
constraint translated to initial conditions at some high energy scale. 



5.3 Regime 3: 

Here \k\ <C |A|, g and we consider as before two regimes for s. However the discussion will 
be much more involved. We give hereafter the main steps. 



(a) s < m, A: 

Since \k\s <C m, A in this case, it is reliable to expand in the small parameter 
ns/M so f t . Adding up Eq. ( |5.1|) divided by h\ to Eq. (|5^ ) divided by h 2 and using 
Eqs. (|5.3| , |5.5|) to eliminate the combinations hih 2 and h\ + h 2 on finds 



with 



KS 



2^- x f L - r] + O((-—Y)=0 



ill 



oft 



X 



\k M\A\ 



A K 



g 2 m 2 s + §M 2 z K A x m 2 + §M 2 ' 



A 2 



1)M 



m 



(5.19) 



(5.20) 



Thus 



fi 



± 



]jl + 0(x) 



provided rj > 0. 

On the other hand, Eq. ( |5.9| ) yields 



(5.21; 



%s 2f3 Ml) = 
9 



Consistency between equations ( |5.22|) and ( p.!9|) requires 



(5.22) 



1 A 2 K 

A = -^-(2m 2 + 77 + (c 2/3 + — (1 - c w ))M 2 z ) + 0(x, - 



on top of Eq.( |5.21 ). Furthermore, combining Eqs.(|5] 
EWSB relation 



(5.23) 

to retrieve the familiar 



M 2 Z _ [m{ + /2 2 ) - (m 2 + fi 2 )tl 



(5.24) 
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and using Eq. ([5.21|) one gets a further correlation between mi,m2,A\ and m s in 
the form 



where 



rrin 



\ 2 

g 2 



l)Mf 



X)c 2 p 



X 



mi 



(5.25) 



(5.26) 



Finally, equating the values fi 2 from Eqs. ( |5.21 , |5.23|) and using Eq. (|5.25|) to eliminate 
rri2 one finds the necessary correlation among the mass parameters mi,m s , A\ and 
M z , for given t p {= j*), 



X 



'1± 



C2/3 



(5.27) 
(5.28) 



where 



Xl = l(2ml + {^{1 

tp g 2 



C2/3) + c 2f5 )M 2 z ) 



(5.29) 



Equation ( f5.27| ) is easily obtained by eliminating ft between Eqs.( |5.21"D and ( |5.23| ) 
and solving for the resulting quadratic equation in X. (For simplicity, we skip from 
now on the explicit indication that all the relations are valid only to zero th order in 
small parameters like k/X, n/g, Ks/M so f t .) 

Incidentally, we note here that 77 > as required by Eq. ( |5.21[ ), is automatically im- 
plied by the correlations Eqs. (|5.25| , |5.27|) . Furthermore, m 2 > implies immediately 



< < 1 



(5.30) 



from ( |5.26|) and (|5.27|) . For definiteness we stick hereafter to the phenomenologically 
likely case tp > 1. Then, it is easy to see from Eq. (|5.27| ) that 1Z + < 1 is excluded 
by tp > 1, while 1Z~ is acceptable and leads through Eq.( f5.30| ) to the constraint 



A 2 > 



6 2/3 



(5.31; 



which overpowers Eq.( 5.28 



One can now determine uniquely \x from Eqs. ( |5.23| , |5.20| , |5.25| , |5.27|) : 
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Now as far as X\ ( Eq.( |5.29| )) remains positive at the electroweak scale (which is the 
case if tan f3 does not become extremely large leading to substantially negative m\ 
at this same scale), equation (|5.32j ) shows that \fx\ is a decreasing function of \A\\. 
Then using Eq.( |5.31| ) one gets the following upper bound on fx 2 



£ 2 < f (5.33) 

This constraint leads to an upper bound on the physical lightest chargino mass, 
which may or may not be consistent with the experimental lower bounds. To go 
further in this issue, let us include first the requirement that the electroweak sym- 
metry breaking extremum of the scalar potential should indeed be a minimum lower 
than the SU{2) L x U(l)y symmetric one at vanishing scalar fields, i.e. that 



Vm.iT), ^ 



(5.34) 



Evaluating Eq. (|5.13 ) at k ~ and \x as given by Eq. (|5.32|) , the above condition 
leads to 



4*5(1 + t 2 p)(x 2 - (1 + tftxJAl + (x 2 - (1 + 4) 2 xi) 2 > (5.35) 
where x% is as given before, (Eq. (|5.29| ), 

,A 2 



*2 = + ^)Mf 



r 



(5.36) 



and Eq.( [5.3l|) has been implicitly used. As can be seen from Eq. (|5.35| ), this equation 
leads to a constraint only if x 2 — (1 + < 0. We need to consider two cases: 



i) %2 — (1 + tp)xi > 0: in this case Eq.( |5.34| ) is trivially verified, but Eq.( |5.33| ) 
together with x 2 — (1 + tV)xi > and ( |5.36| ) lead to 



tf<^{- 2 sl p + c 2 w )Ml (5.37) 



Since the chargino sector is identical to that of the MSSM one can study directly 
the effect of the above bound on the mass of the lightest chargino denoted hereafter 
by M +. Using a rigorous upper bound on M + 
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\ 



fi 2 + 



2M 2 , 
1 + tl 



(5.3* 



in conjunction with Eg. Q5.37| ) one gets immediately a tan/3 dependent (and fi in- 
dependent) upper bound on M +. Confronting this upper bound with the present 

Al 

experimental lower bounds from the LEPII exclusion analyses |^2, 23] one finds typ- 
ically that for tan j3 > 1 only a small range of tan f3 values close to 1 are possible. 
For example, taking A = g one has 



M + < 



1 + t 



2 M^ 



+ 



1 + 

At 2 



Mf 



(5.39) 



Comparing for instance with the experimental analysis of |23[ (Fig 5 therein) one 
finds that only a very small range, 1 < tan/5 < 1.3 — 1.6, is consistent with our 
regime in this case. Of course one should keep in mind the model dependence of 
the experimental analyses. However we should stress that since in the regime under 
consideration A is relatively large (~ g), we are in a configuration which is safe from 
significant unconventional signals due to displaced vertices that can occur in the 
(M+1)SSM |J. The comparison with MSSM-based experimental analyses is thus 
fully consistent here. 



ii) X2 — (1 + tg)x\ < 0: In this case the stability of the electroweak symmetry 
breaking minimum requires that A\ be bounded from above at the electroweak 
scale. This upper bound A\ can be easily read from Eq. fl5.35D l 7 ]. However, A\ < A\ 
will lead, through Eqs.( |5.27] , |5.26| , |5.29| , |5.36| ) to an upper bound on m 2 s (remember 
that Eq.( |5.27| ) is valid for 1Z~ , see the discussion following Eq.( |5.30| )). Working out 
the algebra one finds 



ml < 



A 2 



Mf 



g 2 4m 2 



1)M 



1 + ^) 2 



((*; 



9 



(5.40) 



This inequality shows an anti-correlation between m 2 and m\ which suggests that 
universality assumptions between the singlet and doublet soft scalar masses may be 
disfavoured. To see this for any value of tp one should plug the running expressions 
for m\ , m 2 given in Eqs.( flB.18j ), flB.20|) ) and evaluate the auxiliary functions given in 
Eqs.( (|3.5|) - ( ^.9|) ) and in Appendix B (approximating them for instance to their first 
order iteration). Instead, let us give here a simpler and more qualitative argument. 
It is clear from Eqs.( |5.33"| , 5.37]) and from the dependence of x\ on tp (Eq. fl5.29p ) 
that the higher the experimental exclusion bound on M v + the less favoured are the 



7 which in turn leads to a phenomenologically harmless lower bound on \jj,\, \p,\ > \p,(A\)\, rather than 
an improved upper bound like in i). 
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large tp configurations. This conclusion is reinforced by the fact that for a given 
initial condition m§, the running m\ decreases faster for larger tp, thus driving x\ 
to small values. So let us concentrate only on the small tp region. In this case the 
running m 2 s and m\ are well approximated by the simple analytical solutions given 
TTJ, namely m\ ~ m\ + M^j2,m 2 s ~ m\, in the vicinity of the infrared effective 



m 



fixed point for small tan/3 (~ 1) and assuming universal initial conditions for the 
scalar and fermion soft masses (denoted respectively by mo and Mo). Using these 



relations in Eq.(5.40) one finds 



ml < + + Mf) + 0(tp - 1) < M + (tp - 1) (5.41) 

that is m\ is forced to be rather small (eventually vanishing) for large M5 2 , the latter 
being required by experimental lower bounds on M x + . A way to avoid such an anti- 
correlation is to relax the universality between the singlet and doublet Higgs soft 
masses. 

(b) s > m, A: 

In this case, no generic statement can be made about the size of ks. If the magnitude 
of A is such that M so f t /(Xs) is small then Eqs. (|5.6| , |5\7| ) can be cast in the form 



hi = ||[l + ^^ + 0((^) 4 )] (5.42) 



hi = ^[l + ^^) + 0((^f)} (5.43) 



2^ L 2tf ' " ft 

Feeding the above equations back in Eq.(p.l|) one gets 



A 2 + ? + 5 + 0£(^)\ (^r) x M = (5.44) 
jl A // // 2g z 



where 



7 = -A x {m\-m\) (5.45) 

S = ^-M 2 z + ml + ^-(m\-ml) (5.46) 
T 2A 

In the deep regime 3, |ac/A| 1, thus 5 is positive (except for very large tp where 
m\ can become negative at the electroweak scale). It then follows from Eq.( p.44[ ) 
that this regime can not be dynamically consistent in the case at hand, i.e. as far 
as A is not too small in this regime so that s ^> m, A implies very large |A|s(= jl). 
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Finally let us stress that Eq. (|5.44 ) is more general than in the context of regime 3 
(the condition k <C g was not used in establishing this equation). Furthermore, \8\ 
is by definition of order Mg t t whatever its sign. Thus if s happens to be extremely 
large compared to the soft masses, then Eq. (|5.44|) forbids ft to be very large too, i.e. 
A cannot be large (~ g say). We thus have a further ingredient in understanding 
purely numerical studies were indeed very large s required very small A (see for 
instance table 1 of reference IITll). 



5.4 Summary for regime 3 

Regime 3 is the trickiest. We have shown that an s much lager than the soft parameters 
is forbidden unless the hierarchy \k\ <C |A| < M so ft/ s is realized. On the other hand, in 
the case where s is of the order of the soft parameters the stability of the EWSB vacuum 
has to be invoked. In configurations where this stability is automatic, fx leads to light 
charginos inconsistent with the present experimental limits (except for a small window 
1 < tan (3 < 1.3 — 1.6 which will also be closed by a few GeV improvement of these limits). 
When the stability of the EWSB vacuum is not automatic, then the generic price to pay is 
small tan (3 values and either a relaxation of universality between the singlet and doublet 
Higgs soft scalar masses, or a small universal soft scalar mass mo anti-correlated with a 
large universal soft gaugino mass Mq to ensure consistency with experimental limits on 
the lightest chargino. Clearly future improvement of these mass limits will reinforce the 
above conclusions. 



6 Conclusion 

We have made an extended analytical study of the scale evolution of the various basic 
parameters, as well as of the spontaneous electroweak symmetry breaking and the dy- 
namical determination of the \x parameter in the (M+1)SSM. In particular, we identified 
four regimes of effective infrared fixed points behaviour corresponding to various relative 
magnitudes of the two couplings A and k that enter the gauge singlet scalar sector. These 
regimes correspond to the configurations of least sensitivity to the initial (GUT-scale) 
conditions of most of the parameters. We have determined analytically this sensitivity 
and shown how it generalizes the MSSM case. Furthermore, the analysis of electroweak 
symmetry breaking (which did not require any numerical scan over the parameter space) 
showed that some of these regimes, taken in a wider sense, are generically excluded by 
negative experimental searches for charginos or by purely dynamical considerations, and 
that the others lead to very large singlet vacuum expectation values which can be reduced, 
however, if some amount of non-universality of the soft parameters is allowed. 

The general analytical solutions given in this paper do not rely on any model-dependent 
GUT-scale assumptions. In practice they lead to approximate solutions to the RGE's 
(with controllable convergence), in analytically closed forms and for any value of the tan /3 
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parameter. They thus allow a precise study of the dynamical properties of the (M+1)SSM. 
Besides, they are readily generalizable to further extensions of the Higgs sector. This pro- 
vides a basis to study the genericity of these properties beyond the (M+1)SSM and to 
gain a more thorough understanding of the sensitivity of the Higgs sector phenomenology 
to specific underlying supersymmetric models. 
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Appendix A: The coefficients a^, b x and c& 



In this appendix we define the coefficients and introduced in the RGE (|3TT 



<ht = 6 


Ofet = 1 


a Tt = 


a\t = 3 


a Kt = 


a t b = 1 


abb = 6 


a T b = 3 


a\b = 3 


a K b = 


a tT = 


a bT = 1 


a TT = 4 


axr = 1 


a K b = 


a t x = 1 


ab\ = 1 


a T \ = 1 


axx = 4 


a K x = 6 


a tK = 


dbK = 


a TK = 


axn = 2 


a KK = 6 



bi = 11 



n — 13 


Cfti - | 


C T 1 — 3 


CXI = 1 


c«i = 


Ct2 = 3 


C&2 = 3 


C T 2 = 3 


CA2 = 3 


c k2 = 


r - 16 

Q 3 - y 




C r3 = 


CA3 = 


c k3 = 



Appendix B: Exact solutions for the A's and the El's 

Hereafter we give the (recursive)-equations defining the auxiliary functions which enter 
the A's and E's. 



The A's 



c t 



E t drj 



A° b J u b- Ju b e b , yoAlfux - Juxex 



6Y b °ju b 



:ba) 
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1 dE b v0 A°fu t - Ju t e t A° T Ju T — J u T e T 



+ Y } 



E h dr) 1 l + 6Y t °ju t T l + 4Y r °jX 
o A \J u \ ~ J u \ e \ 



l + 4Y°fu x 



1 dEr _|_ 3 yO ^°/ M fc ~ /^ e ft + yO ^A / M A - / U X e X 



E T dr) b l + 6Y b °ju b l + 4Y°fu x 



1 Ag/^ - Jw f e f Ag/n 6 - Ju b e b 

6X " £ A dij ' l + GYffut +6b l + 6Y°fu b 

_ v q ^a/ m A / ^A e A 

6k "^ ^ 1 + 41^/« A 

where the variations of E k should be taken at r\ = fj = and are given by 

1 dE k 



The E's 



+ 



1 d 2 £ t 2 d-E* 
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0\2 



+ 7(y 6 °) 



i + 4y?/« A 

A b Ju b - Juhet, 1 ' 



l + 6K?/u 6 



0\2 



+ 5(F A U ) 



^a / u x - I u x e x 
' l + 4Y°fu x 



n < -o I u b~ I u b e b ) (A x Ju x -J u x e x ) 



(l + 6n°/« 6 ) (1+4F a °/«a) 



1 d 2 £ 6 + 2 d£ fe 
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l + 6Y t °Ju t T l+4Y°fu T 
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+ Y, 



A^Jjuxj-Juxex 
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- Y 



l + 6Y t °Ju t 
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l + 6Y b °Ju b 



A x I u x - J u x e x 

i + 4r A °/ MA 



(g + (A°) 2 ) J^ - 2A° J Mfe e b + J ^ 
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1 d 2 £ A 2 d£> 
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27(Y t 
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(B.10) 



= 1 d 2 E K | 2 dE K 
E K dr]drj E K drj 
(X°x + ( A °x) 2 )Jux 



+ 



6 



60(F A ) 2 



l + 4Y°fu x 
-2A x Ju x ex + IuxCx 



A \Ju\ - Ju x ex 
l + 4Y x °Ju x 



(B.ll) 



1 d 2 E k 
Ek drjdr] 



v =0,fj=0 



t 2 C£c kt a t M?f + 2tY.c ki a l {M») 2 - ^^(Mf) 2 (B.12) 



8=1 



i=l 



i=l 



The m's 

Due to linear relations |TB| which follow from the RG equations ( |3.1| ), ( ^.25| ) (dropping 
out momentarily the trace term S in the latter equation), one can express the soft masses 
for squarks, sleptons and Higgses in terms of the in the form 
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, 2 , 27/i-81/ 2 



< = Wr+ J 96 ' (R13) 

45(E f - g) + 3(E fe - g) + 6 (E T - S°) - 27(E A - £°) + 9(E K - S°) 



96 



o , 2 , 21/i - 63/ 2 



< = Kb) 2 + Qg (B.14) 

3(E f - g) + 45(E b - g) + 6 (E r - E°) - 21(E A - E° A ) + 7(E K - E°) 



96 



- 2 Q = K) 2 + " 24/l 96 + 72/2 (B.15) 
24(E f - g) + 24(E b - Eg) - 24(E A - E°) + 8(E K - E°) 



96 



, n . • . — 33/i + 99/ 2 , n .. ,, , 

m L = {ml) + — (B.16) 

9(E f - g) - 9(E fc - Eg) + 30(E T - g) - 15(E A - E° A ) + 5(E K - g) 

96 



mi = + (B, 7 ) 

18(E f - g) - 18(E fe - g) + 60(E T - g) - 30(E A - g) + 10(E K - g) 

96 



To! = (ml) + — (B.18) 

-27(E f - g) + 27(E b - Eg) + 6(E T - S°) + 45(E A - £°) - 15(E K - E°) 
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> , n , > 3/i + 9/2 . . . , n 

™2 = K) + gg (B.19) 

27(E f - g) - 27(E b - g) - 6(E r - g) + 51(E A - g) - 17(E K - g) 

96 

™ 2 = K) 2 + ^(s K -s°j (B.20) 

where 



6j V (l + Moi) 2 / 
To obtain the complete solutions one should add to each of the above equations the 
corresponding trace term contribution, see Eqs. (|3.29|) , ( |3.30|) and the discussion in section 
3. 
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Appendix C: The Y's in the IRQFP regimes 

We give in this appendix some details about the derivation of Eqs.( |3.14l) - (|3.16|) . For 
later use we define 

Y° 

<oo (C.l) 



when it makes sense, namely when Y t ° and Y® go simultaneously to infinity with a fixed 
ratio rj. 

We will show inductively, in the same spirit as JIfi|], that in such a regime Ui — > uf with 



FP 



(Y t ° 



(C.2) 



where the p^s are fixed numbers which we will explicitly determine, and the uf p, s are 
initial conditions independent and will be defined implicitly through equations ( |C6| ) - 

To proceed we consider finite order iteration approximations to Eqs.( |3.5| )- (|3.9I). To 
obtain the (n + l) th order approximation to Mj, denoted w 4 - n+1 ' ) , in terms of the u\ n \ one 
makes the formal substitutions Ui — > w,- n+1 ^ and — > respectively on the lefthand 
and righthand sides of Eqs.( |3.5|) - (|3.9| ). If for any given one of the four IRQFP regimes 
of section 3.2, the u^s have the behaviour of Eq.( |U.2j ) at the n th order, i.e. 

,FP(") 



U 



in) 



U 



(n) 



(C.3) 



where the scale dependent functions u FP ^ are Y t ° independent but possibly r b ,r T ,r x ,r K 
dependent, then the same is true at the (n + l) th order with the following recursive 
equation for the p^s 
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(C.4) 



where we define 



A ij [q,x} = {S qi + 5 qj )e[x] (C.5) 

the <5's are Kronecker's and 9 is the Heaviside function. Equation ( |C4| ) describes com- 
pactly the four IRQFP regimes labeled by q = 1, 2, 3, 4. The 9 function is here to account 
in general for the fact that some of the p^s can become larger than one at some iteration 
order. It is important to keep this point under control since in a regime where some is 
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Table 1: Powers entering Eq. ( |C.2j ) in the various IRQFP regimes 

very large, one hasQ 1 + Y t ° J* ~ Y? /„* if pj n) < 1 and 1 + Y? J* u[ n) ~ 1 if pS n) > 1 
( the critical value pf 1 ^ = 1 is never met). Thus the proof of convergence to a unique form 
for the w°°'s ( i.e. definite limiting values for the p^s and for the functions uf p ) needs some 



care. Direct numerical inspection of the iterations in ( C.4|) shows that only p^ can be 



either bigger or smaller than one. For regimes 1,2,3, p^ is found to be constantly smaller 
(or bigger) than one at each iteration. Regime 4 leads to some slight complications: p^ 
and p^ take alternating magnitudes, > l,p^ > implying < l,pj-. n+1 ) = 

and vice versa up to the first six iterations, but then stabilize with p A > 1, p^ = for 
all n > 7. Armed with this numerical information, one can then consider Eq.( |C.4|) for 
any n > 7 to prove recursively that for all four regimes < pf^ T K < 1, while < p^ < 1 
for regimes 1,3 and p^ > 1 for regimes 2, 4. Eliminating the 9 functions correspondingly 
and solving Eq.( |C.4| ) for pf 1 = p$ — Pi, we find the limiting values given in the Table 
1. 

The same result can be derived through the more systematic method of solving directly 
the linear system (|C.4|) for pf 1 ^ = pf^ = Pi, keeping formally the 6 functions and not 
relying on any prior numerical information about the iterations. Then, checking the 
consistency of the solutions for all possible combinations of the p^s being greater or 
smaller than 1, one delineates the acceptable ones. It turns out that there is only one 
consistent solution for each regime leading to the values given in the table. 
Knowing the p^s one can now extract the itf p 's, 



1*1 regimes 1, 2 and 4 
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^provided that t is also large enough. 
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-t— J r regimes 1, 2 and 4 

Et 

(6r fc )i(4r A )I(/ « FP )?(/« FP )3 
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— t 1 1 1 1 t regime 1 
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Si 
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(6)3(6r t )3(4»y)*(/«* I, ) !! (/ « FP )?(J < P ) 


(6)i(6r 6 )^(4r T )M6r K )i(J < p )2 (/ (J < p )4 (J < p ) 



regimes 2 and 4 



{1 regimes 1, 2 and 4 
(CIO) 
pp 3 regime 3 

Finally let us note that p K is non vanishing (but small) only in regime 3. Consequently, 
the evolution of Y K will be very slow in this case, and a IRQFP regime will be theoretically 
obtained for very large values of Y" K ° (or very large values of r K ), as it has been seen in the 
numerical analysis of section 4. 

Appendix D: The soft parameters in the IRQFP regimes 

In this appendix we discuss in some detail the sensitivity of the soft couplings and masses 
to their initial conditions in the IRQFP regimes. Keeping in mind that one can derive the 
results directly from the Yukawa couplings discussed in the previous appendix through 
the substitutions Eq. fl3.10D ) it is instructive to make a direct study starting off directly 



from the analytical forms in the soft sector. 



The A's 



To understand the behaviour of the soft trilinear couplings in the various IRQFP regimes 
let us study first that of the auxiliary functions given in (|B.1|) - (|B.5|). To start with, 



we take for illustration the case of e t . Denoting by e°°'s the limit of the e^'s when Y® bT 
tend to infinity, and assuming that to is sufficiently large so that 1 + 6Y b ° / *° Uh is well 
approximated by 6Y b ° J^uf, Equation ( |B.1[ ) reads 



e ~(t ) = -^(t ) + -A° b - ^ "f eT + y o ^Uo°^-/o°<er ( } 

The above equation is valid for any of the four regimes considered in section 3.2. To 
write more specific forms for each regime one uses Eq. ( |3.17| ) together with the correspond- 
ing values of the powers given in the table of Appendix C. One then obtains 
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E t drj Q$uf 6 
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l A °x o 

_l_ _a regime 3 



DC , .Xl- 

A fc A 



(D.2) 



In deriving the above equation for regimes 1, 2 and 4 we have made the implicit 
assumption that the magnitude of remains under control so that / u^e^ — > with 
growing y°. 

Equations ( |B.2| ) - (|B.5| ) yield along the same lines the corresponding e°° in the various 
IRQFP regimes: 
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regimes 1, 2 and 4 



+ regime 3 



(D.6) 



where we used Eq. flC.10D when convenient. At this level we should stress that the particu- 
lar structure of the above equations will actually allow to solve exactly for the dependence 
of the e°°'s on the initial conditions A®. Indeed, on one hand the A^s enter linearly the 
inhomogeneous parts of the integral system of equations ( p.2|) - ( p.6|) with t indepen- 
dent coefficients. On the other hand, all the integrated parts which induce iteratively a 



dependence on the A's are of the form 



so that any substitution therein of the form 



e°° — > J2j el A®, where the c-'s are scale independent constants, will yield again the same 
kind of dependence. It is crucial that all the coefficients multiplying the A's remain scale 
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Table 2: Resummed linear weights a\ of the A°'s in the various e°°'s, i.e X)j Q^j 1 C e°°. 
The unbracketed numbers are common to regimes 1,2,4 and the numbers between brackets 
correspond to regime 3. 

independent at any order of iteration of the integral system of equations. It follows that 
one can re-sum the ensuing infinite series giving the numerical coefficients of the v4°'s. 

Equivalently, one can replace formally everywhere J / J by ef 3 and solve the resulting 

J U i 

linear system in the e°°'s to obtain the exact linear dependence on the A's. The e°°'s take 
then the following form 

er = ef p + 5>M° (D.7) 

3 

where the ef p 's are completely independent of the ^4°'s. The various coefficients a\ 
corresponding to the full resummation have been summarized in Table 2. Note that 
there are in Eq. (|D.5| ), regimes 1, 2, 4, two terms not respecting the convenient structure. 
Nonetheless, due to a conspiracy among the various regimes, these terms do not invalidate 
the procedure described above for the determination of the coefficients aj, even though 

is scale dependent in regimes 1, 2, 4 (see Table 2). 
As for the ef p 's, their defining equations are the same as those for the e?°'s Eqs.( |D.2"|) 
- ( |U.(j| ) with all the A°'s dropped out (the latter cancel automatically when Eq.( [L).7| ) is 
used). 

We have now all the ingredients to determine the dependence of the various A°° on 
their initial conditions. Since in all four regimes Y® b T are large A°l t b T take the form 

f n oo e ao fu Fp e FP 

— e°° + j^J — ~ e i + j u jp , see Eqs. (|3.12|) , (|D.7|) , and are completely insensitive to 

the initial conditions A® similarly to the MSSM case [pf] . In contrast A\, A K will be both 
sensitive to initial conditions in some regimes. Using Tables 1 and 2, equations 



(|C.2[) , (p.7|) one finds in the regimes j L' ^' ^ 
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2Y°t 



^K. [1] 

^°=^f-<! /j 1 ^ A° b -\ * + ^ ^ A° x -\ Ul [2,4] (D.8) 




[3] 



and 



A™ = Al p + { % A° + { I A* + { % A* - { I A° + ( +6¥k °[2,1] ^ (D.9) 

I 31 I 31 I 31 I 31 \^ j^O [3j 

where Af F is denned in Eq.( |3.21| ) 

Hence, when one goes beyond the MSSM it is important to distinguish between e FP 
and e°°, the former being useful intermediate functions to define the initial condition blind 
parts of the Ak in some IRQFP regimes. 

Finally let us note that the sensitivity to the initial conditions in Eqs.( p.9|) , ( p.8|) does 
not imply that the physics is no more blind to these conditions. As was already stressed 
in section 3.2, the Yukawa couplings that multiply A\ iK in the Lagrangian are vanishing 
in the corresponding IRQFP regimes so that the expected screening properties are always 
recovered at the level of the physical quantities. 



The E's 

The study of the auxiliary functions £j is technically more complicated than that of the 
ej's. The discussion goes along the same lines as in the previous section in that the scale 
independent initial conditions contributions can be easily resummed, but there also appear 
scale dependent contributions in the £j's which we should discuss carefully. To illustrate 
the case let us consider for simplicity the top/bottom sector switching off momentarily 
all other contributions. In this case Eq. (|B.7|) reads in the limit Y t ° b — > oo 



+ 



1 d 2 E t 2 dE, 
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and similarly for Eq.( [B.8| ), 
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(D.10) 



6° = £° 



(d.ii; 



Apart from the constant terms (E" t + (A bt ) ) which can be fully resummed there is a 
non trivial dependence on the initial conditions in the above equations, either explicitly 
or implicitly through e^t, which would be difficult to handle in the iterated system (|D.10| ), 
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(|D.11| ). Fortunately the corresponding terms will actually all cancel out. To see this, it 
is convenient to make first the following change of variables in Eqs.( p.lO|) , ( p.ll|) : 



C = &°° + P k ^l + pM) 2 + 2p k 3 A° k e? (D.12) 

(k = t,b), where p\-,p\ are arbitrary constants which can be taken equal to 1. In the 
ensuing equations for £ t , £ b we use Eq.flD.71) to extract all the dependence on initial 
conditions in e^. In the resulting equations we substitute for in £ 4 and vice versa. 
If now pg is chosen in Eq. ( p.l2j ) such that p\ A^ b = Y^jO^bA^, see Eq.( p.7| ), and using 
Eqs.( p.2|) , (p.3| ) for the reduced top/bottom system, one finds that all terms that contain 
initial conditions and are scale dependent cancel out from the defining equations for £ t 
and 5 []. 

Thus we are lead to a situation to the one for e t , tb of the previous subsection. One can 
then similarly define 



= e k p + e 7^? + e iu m (d. 13) 



'■j 



and determine the 7's by solving a linear system of equations. All the ingredients described 
here generalize to the complete t, b, r, A, k system albeit tedious algebraic manipulations. 
We will not write in this appendix the rather complicated dependence of the £°° on the 
GUT scale values of the soft terms. We give here directly the dependence in the Si's. 
A useful remark is in order to understand this dependence. A direct calculation shows 
that the expression of S^ t6 , Eq.( |3.12| ) in the limit Y® — > 00, is invariant under the 
substitutions e°° — > = e°° — p\A®,t^° — ► £i,£j — ► £f p where £i,£f p are as defined 
in Eqs.( p.l2|) , ( p,13| ). Equating further the free parameter p l 3 to Y,jOi\A^ so that £f p is 



guaranteed to be A , S° independent and ei = ef p , all dependence on the initial conditions 
of the soft parameters drop out from S?°. Defining 

r FPt-FP 

Sf p = (AD 2 + 2ef p A pp + C FP - (D-14) 



one has 



and 



S^ Ar = S pp (D.15) 
S~ = S pp + /3iS; (D.16) 

S^Sf+^S? (D.17) 



9 Thc specific values of the ct\ b (j = t, b) are of course crucial for this cancellation. It should be clear 
however that in the simplified example of the top/bottom sector we consider, these values differ from 
those in table 2. 
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where the E£ p r 's are initial conditions independent for all four IRQFP regimes while 
E^'s can obviously still have some dependence in the regimes where F A ° or Y® are not 

' [1,2,4] 



large. One finds for the latter in the regimes 



[3] 



A 



s FP - 





f 4 [i] 

§E° [2,4] (D.18) 
[3] 



and 




dl [1] 
[2,4] 

K [3] 



(D.19) 



with 



A 1 + QYH 



E„ — 



1 + QYH 



(D.20) 



and 



The C F read 



0\2 



(1 + 6F K °t) 2 1 + 6Y°t 



(D.21) 



ef P = 



1 d 2 £t 



E t dridfj E t dr] 3 J « 



pp 
b 



+ 
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(D.22) 



?6 



1 d 2 g b 



1 rfgj, 
B 6 djy 



1 7 f .K P *r P 

+ 36 ^ J u p 



tFP _ 1 d 2 g T 
E T dr]drj 



T 16 



f« FP e F1 
2 

+ * 



,K p ff p 



12 



1 dC T J u I Pe I P 



E T dr] J u 



FP 
b 



2/<P 



+ 



FP 

6 



4 J u FP 



(D.23) 
(D.24) 



1 d 2 E x 
Ex drjdJr] 



J_dE± 

Ex dr. 



_/< P « 
2 f u , 

/»r p er p N 2 



FP 
_j 

FP 

t 



/< P 



+ 



FP.FP 
6 e i> 



FP 

6 



L< p «r 



FP 
6 



/ 

,K p e 



fK Fp e FP 



f< p *r 



FP 

'6 



+ 



l/«f p ep , /uPef p 



6 °b 
FP 



+ 4-7-, 



/< p e FI 



16 



f» Fp e FP 



+ 



l/< P < P /» FP ^ P 



(D.25) 
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FP 



E K dr]dr] 

It turns out that in regimes 1, 2, and 4 one has df = 



(D.26) 



while in regime 3, we have 



2 dEtfu 



FP„FP 
A e A 



E t dr] AJu 



FP 
A 



FPfFP 
A 



AJu 



FP 
A 



+ 16 



+ 



12 /< p 



(D.27) 



2dE b Jufef Ju^i 



,FP c-FP 
^A 



FP 

A 



AJu 



FP 

A 



16 V /u 



FP 
A 



1 /« 



FP^FP 
t 



FP^FP" 



12 JX 



FP 



+ 



8 /< 
(D.28) 



FP 



FP^FP 



2 dE T Ju{ F e 



x 



E T d V AJuf 



i FP 16 



FP 
A 



+ 



1J> 



FP^FP 



FP„FP 
A e A 



4 /« 



FP 

6 



FP 

A 



(D.29) 



For d*x we find 
• regime 1 



dl 



! <^'„ y<) ./ 9y»J_ii I -icfycm / J e « p V 

E K dr, 1 k l+6Y°t l+6Y°t xu ^k/ ^ l+6Y°t y 



r fp 
4 dE K -^0 J e « 

« 1+6Y 



« l+6Y°i 



PP„FP 



FP 
t 



+ 2 



r r . pp 



FP 
b 



(D.30) 



• regimes 2 and 4 



6t 



2 ,IL K Je^ 



E K dr] 



3t 



.k p /v- 



FP 

t 



+ 2^ 



3t 



FP FP 



£ p 



+ 



t 

/< P e FI 



(D.31) 



• regime 3 



(D.32) 



and for dt 



<2r? J'u 



3 d£ K J « FPe A P 



FPtFP 
3 e A jA 



FP 
A 



7^ 



FP 

A 



FP 



regimes 1,2, and 4 



regime 3 



(D.33) 
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The m's. 



When we develop the expressions (B16-B23) in the different FP regimes and replace the 
S t by their dependence on the soft masses mj, we find solutions of the form 

m \ = ( m i)m» + ( m ^M» + (rn^Tr (D.34) 

where (mf)^ is given by (|3.30|) and the (mf) m o are written below. But, there is no 
compelling exact analytical solutions for the {mf) M o. For the regimes 1,2 and 4 : 

, ,o 85 i)2 17 q 2 10 o 2 5 n 2 5 q 2 15 q2 17 n2 

(mn) m o = m n> m,r m n , H m T H m Fo m, m 9 

v vv™ 122 q 3 12 2 3 61 3 122 L 122 £3 122 1 122 2 

/ \2 17 n 2 40 n2 4 n 2 1 n2 1 n 2 3 n 2 21 n 2 

m-f L = m n H mf> H "n^n wip H ?nY ^ 9 

/ \2 20 n 2 ^ n 2 37 n 2 6 q 2 6 n 2 18 n 2 

(m B ) m0 = -^j-m Q3 + ^-m^ + — m D3 + — m L + — m B3 - — m x + 

. .n 15 n 2 3 q 2 9 n2 87 n2 35 q2 17 n2 3 n2 

(mr) m o = m n , m (/ , H m n „ H m r m Po m, m 9 

v Lj m 122 q 3 122 r/ 3 t 61 3 122 L 122 Ei 122 1 122 2 

. n 15 n 2 3 Q 2 18 |i2 35 n2 26 n2 17 n 2 3 q2 
mpLo = — mn,, m TJri H m n . H ti F) rn? 

V *>Jm gl Q 3 gl C/ 3 ^ 61 A> gl i ^ 61 ^3 61 1 61 2 

/ \2 45 n 2 9 q2 27 o 2 17 n2 17 q 2 
mi Lo = m ni H m lu m n , m T m Fo 

V 17m" 122 Os ^ 122 U 3 gl D 3 122 L 122 £ 3 



4 2 

— m 9 
61 2 



, ,o 51 n 2 63 n 2 6 n 2 3 n 2 3 n 2 9 n 2 59 (|2 

m 9 in = Tin "mli H mi f^w H tiV H m 9 

V 2J m o 122 Q 3 122 C/ 3 T- 61 D 3 122 L 122 £ 3 122 1 



51 Q 2 9 Q 

m, H m 9 

122 1 122 2 

,7l 9 

122 - 



!m% 2 +396Y0m% 2 -132Y0(A K ) 2 . 1 

; (i+39W) 2 re S ime 1 (D.35) 

m° s regime 2, 4 

In the regime 3. we obtained 



o 39 q2 11 g2 6 q2 5 q2 5 q2 5 g2 17 (l 



TioL,o = — m n mjr m n H m T H m w mV m 9 H 

, >o 11 n 2 19 rj 2 1 n2 X n 2 1 ri2 3 r\2 9 r\2 3 ri2 

mT m o = m n H m TT H m n m T m p H m, m 9 H m q 

\ Dm 31 & T 31 u 3 -r 31 d 3 31 31 B3 31 1 31 2 31 5 

/ \2 12 n 2 1 n 2 l^ n 2 8 n 2 8 n 2 8 n2 10 n 2 7 n 2 

me Lo = m n „ H m r r H m n , H m T H m p , H m 9 H mo 

. , 9 5 n 2 3 q 2 4 q 2 131 q 2 55 q 2 7 q 2 1 q2 5 q 2 

{m-Ll^o = — m n „ mjr, H m no H m r m Pi m, H m 9 H 

V ^Jm g 2 Q 3 62 % -r 31 i? 3 ^ lg6 i 18 g ^3 g 2 1 lg6 2 -r g g S 
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. x2 5 q 2 3n2 8q2 55 n 2 38 n 2 7 n 2 1 n 2 10 o 2 

jtie Lo = — Tn n Tnfr H m n m T H m p 777.1 m 9 H 

£jm 31 Q3 31 Us 31 Ai 93 L 93 Ba 31 1 93 2 93 5 



, ,n 15 ri 2 9 r.2 12 q 2 7 q2 7 r.2 21 n2 1 n2 5 n 2 

(rni „o = m n H m,, mn m> mVr H mV m, m q 

v w 62 q 3 -r g2 c/a 31 d 3 g2 l 62 Ss 62 1 62 2 31 s 

i \l 1' n 2 27 n 2 5 n 2 1 n 2 1 n 2 1 n 2 71 r.2 17 

m 2 L,o = m n rriir H rrin H m> H mV H m, 77 

, ,n 16 ||2 9 r. 2 7 Q 2 10 r> 2 10 n 2 10 r.2 34 r.2 32 n 2 

mcLo = — m n H nxij H mv, H m T H mV m, H 

v A; m o 31 Q 3 ^ 31 % -r 31 % ^ 93 93 ^ 31 i 93 2 93 5 

(D. 
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